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Abstract 

In quantum electrodynamics, photon-photon scattering can be the result of the ex- 
change of virtual electron-positron pairs. Effectively, this gives rise to self-interaction 
terms in Maxwell's equations, similar to the nonlinearities due to polarization in non- 
linear optics. These self-interaction terms vanish in the limit of parallel propagating 
waves. However if the modes generated in bounded regions are used, there will be 
a non-zero total effect. We show that stationary two-dimensional light bullets can 
form in guiding structures, due to the balancing effect of quantum electrodynamical 
vacuum nonlinearities on dispersion and diffraction. These light bullets are unstable 
and exhibit the possibility of self-focusing collapse. The consequences of our results 
are also discussed. 

PACS numbers: 12.20.Ds, 42.50.Vk, 42.65.T 



1 Introduction 



According to quantum electrodynamics (QED), the non-classical phenomenon 
of photon-photon scattering can take place due to the exchange of virtual 
electron-positron pairs. This is a second order effect (in terms of the fine 
structure constant a = /Aneohc ~ 1/137), which in standard notation can 
be formulated in terms of the Euler-Heisenberg Lagrangian density [1,2] 

=S^EH = £o^ + e(4^' + 7^'), (1) 

where ^ = 2a'^eln^ / A^mlc> , ^ = \{E'^ - c^E^), ^ = cE-B, eis the electron 
charge, c the velocity of light, 27ih the Planck constant and rrie the electron 
mass. The latter terms in expression (1) represent the effects of vacuum po- 
larization and magnetization. We note that ^ = 5f = in the limit of parallel 
propagating waves. It is therefore necessary to use other waves in order to 
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obtain an effect from tlie QED corrections. Furthermore, it has recently been 
shown that QED vacuum nonlinearities can lead to self-focusing of laser beams 
[3], and that such effects might be experimentally viable within a decade. In 
Ref. [4] it was found that guiding structures can be useful for the purpose of 
detecting QED effects, and a suggestion was presented for an experimental 
setup that gives a measurable signal. 

In the present work we thus consider waves guided by two parallel conducting 
planes, in order to study the possibility of using nonlinear QED effects to 
balance diffraction and dispersion in directions parallel to the planes. We find 
that our system is governed by a 1 + 2 dimensional nonlinear Schrodinger 
equation (NLSE). Using approximate variational methods, which have proven 
to be useful in similar situations (see e.g. [5,6]), we find that the nonlinearities 
can counteract diffraction and dispersion to form 2-dimensional light bullet 
solutions, i.e., solitary solutions to the NLSE which preserve the envelopes of 
the fields. The formation of light bullets was in other contexts (optical media) 
first studied by Silberberg [7]. 



2 Light bullet solutions 



In a medium with polarization P and magnetization M, the general wave 
equations for E and B are [4] 



1 a^E 



V'E = -iio 



^ + ,.V(V.P) + |(VxM) 



(2) 



and 



1 d^B 



V'B 



V X (V X M) + ^(V X P) 



(3) 



Furthermore, the effective polarization and magnetization in vacuum due to 
photon-photon scattering induced by the exchange of virtual electron-positron 
pairs are given by (see, e.g., Ref. [3]) 



and 



P = 2^ [2{E^ - c^S2)E + 7c^(E • B)B' 



M = -2c^^ 2{E^ - c^B^)B + 7(E • B)E 



We consider propagation between two parallel conducting planes with spacing 
xq (i.e., the region < x < Xq is vacuum surrounded by the plates that, as a 
starting point, are assumed to be perfectly conducting). We assume that only 
one TE„o mode {n = 1, 2, ...) is present. To linear order, this gives the fields 



2 



Bj.——ikA sin 



— A cos 

,Xo 




) 



) 



exp[i(A;^ — urt)] + c.c, 



exp[i{kz — cut)] + c.c, 



(4b) 
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exp[i{kz — ut)] + c.c, 



(4c) 



together with ~ cu'^ / — k^ —n^n^ / xl. Here c.c. stands for complex conjugate. 

It should be noted that the beams propagate with a 90 degree angle with 
respect to each other for the selfconsistent solutions obtained in a similar 
situation in Ref. [3]. This would correspond to the relation k = nn / Xq in our 
case. However, we have no such limitation in our derivation. The result of a 
90 degree beam angle of Ref. [3] comes from requiring an explicit symmetry 
between the wave equations for E and B. A similar requirement is not needed 
in our case, since the TE„o-inodes introduced by the conducting planes do not 
have a complete electric-magnetic symmetry as, for example, the electric fields 
vanish close to the planes. What one still might worry about is whether or 
not both Faraday's and Ampere's law, which are equivalent to the two wave 
equations for E and B, are properly solved in our case. However, the equation 
that will be derived for the vector potential is equivalent to Ampere's law, 
which is equivalent to the wave equation for the electric field, and using the 
vector potential we therefore automatically solve Faraday's law. It thus follows 
that our solution for the vector potential corresponds to electromagnetic fields 
that are consistent with the two wave equations for E and B. 

We have expressed the fields in terms of the vector potential amplitude A, 
given by A = (0, A, 0) , where 



using the radiation gauge ((6 = 0). A nonlinear dispersion relation can be de- 
rived by inserting the linear expression for the fields into the right hand side of 
Eq. (3), taking the 2;-component of the resulting equation, and separating into 
orthogonal trigonometric functions. From that equation the coefficients in the 
NLSE can be found. However, here we will present an alternative derivation, 
starting directly from the Euler-Heisenberg Lagrangian (1), which is more el- 
egant and gives the same result. Naturally we must first express in terms 
of the electromagnetic potential, since it is the basis for the original varia- 
tional principle. We take A — A{t, z), and assume A to be weakly modulated 
{\dA/dt\ <S \ujA\, \dA/dz\ <^ |^^|), but we omit the slow y-dependence, since 
such a dependence makes the derivation technically more complicated without 
adding extra understanding. A diffraction term in the y-direction is lost by 
this procedure, but this effect is trivially added afterwards. To lowest order, 
we thus omit the nonlinear terms and the slow derivatives in jSf . Noting that 
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^ = {l/2)[{dA/dty - c^{VAf] and W = 0, we find that the non-oscillating 
part (i.e. the part that does not vanish after integration) is 

^„/.„ . # = h.^ - k^c')\Ar-.n^ ('^) - ^lir^cos^ (^) . (6) 

Z \ Xq / ZXq \ Xq / 

After performing the integration over the region between the plates we find 
that the lowest order part of the action vanishes identically as the dispersion 
relation is here considered to be satisfied. Going to the next order of approxi- 
mation in the Lagrangian we include first order slow derivatives, which yields 
^i/eo = ^o/£o + icu[idA/dt)A* -{dA* /dt)A] -ikc^[{dA/dz)A* - {dA*/dz)A\. 
After variation this leads to an equation where the envelope moves with the 
group velocity. The next order and final approximation includes second order 
slow derivatives and the term in the Lagrangian. After performing the 
,T-integration, dropping (since the action is identically zero because of the 
dispersion relation) and eliminating the second order slow time-derivatives 
using d^/dt'^ ~ v^d'^/dz'^, the final expression for the Lagrangian is 



{^.)^.u:e.[^A*-^A^ 



ikc £o I —A 



dA 



dz 



+ 



4 1^1 ' 
Xn 



(7) 



where () stands for integration over x. Variation with respect to A* leads to 
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£03^0 



^|2^ = 0, (8) 



where Vg is the group velocity and v'g the group dispersion that follows from 
the linear dispersion relation. We have also added the diffraction term in the 
y-direction corresponding to the full amplitude dependence A — A{t,y,z). 
Changing to a system moving with the group velocity while rescaling the 
coordinates and the amplitude according to 



we obtain 



cut/2 
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uy/c 
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da d^a d^a , .n 
which corresponds to the rescaled Lagrangian density 
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where the diffraction term in the ^/-direction is also included. 



From now on we will look for cylindrically symmetric solutions of Eq. (10), 



i.e., a = a{t,p), where 



V 



+ ^2. Equation (10) can then be written 



.da 1 d I da^ 
dr pdp \dp. 



+ \a\ a = 0, 



(12) 



while the Lagrangian density (11) takes the form 



da 



dp 



(13) 



with the action given by ^ = / ^ pdp dr. Equation (12) is a 2-dimensional 
radially symmetric NLSE. Although exact solutions of this equation are not 
available, numerical and approximate analysis give a clear picture of the so- 
lutions. In particular, it has been shown that Eq. (12) allows a stationary 
solution where the diffractive/dispersive spreading of the pulse in the coor- 
dinate p is balanced by the focusing effect created by the nonlinearity. An 
accurate analytical approximation of the dynamics of the pulse-like solutions 
of Eq. (12) can be obtained using direct variational methods involving the 
Lagrangian ^ given by expression (13) and subsequent Rayleigh-Ritz opti- 
mization based on suitably chosen trial functions (see e.g. [6] and references 
therein). A convenient trial function for the present problem is [5] 



aTir,p) 



sech 



P 



fir) 



exp 



mp' 



(14) 



which involves a complex amplitude F{t), a pulse width /(r), and a quadratic 
phase function modeling the phase front curvature. Inserting this ansatz into 
the variational integral and integrating over p, a reduced variational problem 
is obtained for the parameter functions, -F(r), -F*(t), /(r) and &(r). The 
subsequent reduced Euler-Lagrange equations can be rearranged to give F, 
F* and b as explicit functions of the width /(r), which satisfies the equation 



dr^ 



(15) 



where 7 = 4(ln4 + 1)/(27C(3)) ^ 0.29, /(r) = P{t)\F{t)\'' = /2(0)|F(0)|2 = 
Jo, and Ic = (21n2 -|- l)/(41n2 — 1) ~ 1.35J3] Obviously, stationary solutions 
exist when the pulse power satisfies Jq = Ic- Furthermore, the full solution of 
Eq. (15) is 



/(r) = /(0). 1 + 



7 



(16) 



^ Here C{p) is the Riemann zeta function. 
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which shows that the stationary solution is unstable and either collapses to 
zero width in a finite time when Iq > Ic, or diff'racts monotonously towards 
infinite width when Iq < Ic- 



3 Discussion and conclusion 

From Sec. 2, it is clear that the most interesting alternative is Iq > Ic, in which 
case the QED vacuum nonlinearities play a crucial role, leading to a collapse. 
In dimensional units the inequality Iq > Ic roughly leads to Ef^^^k'^rf^^^ > E^^^^, 
where -Einit and rinit are the initial electric field and beam radius, respectively, 
and the critical QED electric field is defined by E^^.^^ = £o/C- Naturally this 
unbounded self-focusing will eventually be saturated by some kind of higher 
order nonlinear mechanism. However, this will not occur before an electric 
field level of order E ~ Ecth is reached, which is the field strength when both 
our perturbative nonlinear calculation scheme and our starting expression, the 
Euler-Heisenberg Lagrangian, breaks down. For such extreme energy densities, 
higher order Feynman diagrams must be included in the QED description, and 
possibly the corresponding physical effects may then counteract the collapse 
scenario, resulting in a saturated beam radius Tgat ~ ^init-E'init/-E'crit- The possi- 
bility to reach such extreme intensity levels is of very much interest. However, 
it is clear that setting up the conditions necessary for vacuum light bullets to 
be formed, is a technological challenge. 

Presently the electric fields that can be supported by the walls before field 
emission takes place are of the order 10^ V/m [8]. For such field strengths we 
note from the above estimates that we must have kri^n ~ 10^^, in which case ei- 
ther the initial beam radius becomes unrealistically large (of the order of km), 
or the wavelength to short (i.e. in the short uv or x-ray regime) for the walls 
to be conducting. Two things should be stressed, however. Firstly it should be 
noted that the experimentally possible field level before field emission takes 
place have increased significantly in recent years [8]. Secondly, in our geom- 
etry there is no normal component of the electric field at the wall surfaces, 
which may allow for a higher (central) field strength than normal, relaxing the 
necessary value on /crinit- A slightly different but related technological ques- 
tion is whether the huge beam powers needed for vacuum self-focusing can be 
reached. This issue was discussed in Ref. [3], where it was concluded that it 
may occur within the next 10-15 years, given the current rate of technological 
improvement. 

A very interesting question from a principal point of view, would be whether 
it is possible to have fully three-dimensional QED-structures which do not 
require any guiding support. However, it seems that this issue cannot be ad- 
dressed within a perturbational approach, and further research is thus re- 
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quired. 



In this paper we have shown that the photon-photon scattering, due to the 
exchange of virtual electron-positron pairs, that effectively gives rise to self- 
interaction terms in Maxwell's equations, can cause optical collapse in vacuum 
if the right conditions are met. More specifically, high intensity electromagnetic 
waves guided by two parallel conducting planes can form light bullets, which 
can collapse if the intensity of the beams are high enough. It is possible that 
higher order corrections to the Euler-Heisenberg Lagrangian can counteract 
such a collapse scenario, thus leading to stable optical vacuum structures. 
Further research is also necessary to find out more about how the optimum 
configurations can be changed, both in the present paper and in previous work 
(e.g. [3]). 
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